Topological invariant in three-dimensional band insulators with disorder 
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Topological insulators in three dimensions are characterized by a Z2-valued topological invariant, 
which consists of a strong index and three weak indices. In the presence of disorder, only the strong 
index survives. This paper studies the topological invariant in disordered three-dimensional system 
by viewing it as a super-cell of an infinite periodic system. As an application of this method we 
show that the strong index becomes non-trivial when strong enough disorder is introduced into a 
trivial insulator with spin-orbit coupling, realizing a strong topological Anderson insulator. We 
also numerically extract the gap range and determine the phase boundaries of this topological 
phase, which fits well with those obtained from self-consistent Born approximation (SCBA) and the 
transport calculations. 
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Time reversal invariant band insulators of non- 
interacting electrons are basically divided into two 
classes: the ordinary insulator and the topological insula- 
tor PHI]. The latter is a novel phase of quantum matter. 
It has an insulating bulk gap and gapless edge or surface 
states. These gapless states are topologically protected 
and are immune to non-magnetic disorder. Recently an- 
other kind of non-trivial quantum phase termed topo- 
logical Anderson insulator (TAI) has been predicted to 
exist in two dimensions (2D) [SHI] and three dimensions 
(3D) [5], which makes the situation more interesting. In 
the TAI phase, remarkably, the topologically protected 
gapless states emerge due to disorder. 

For systems without disorder, the topological phases 
can be characterized by studying the gapless states as 
obtained e.g. from diagonalizing the Hamiltonian in a 
geometry with edges or surfaces. They can also be char- 
acterized by the topological invariants calculated from 
the bulk Hamiltonian, which have been well studied in 
recent literature [10] . However in the disordered sys- 
tems, gapless modes alone cannot unambiguously iden- 
tify the topological phases because they may be local- 
ized in space. Instead, the transport properties are usu- 
ally used to find these extended topologically protected 
modes. Similarly we may also use the topological in- 
variant to characterize the topological phases induced by 
disorder. A question naturally arises: how to calculate 
the topological invariant in the presence of disorder. 

At first glance, it is not obvious how to generalize the 
present methods from translation invariant band insula- 
tors to the disordered systems. Let us first recall the 
integer quantum Hall effect (IQHE) where the general- 
ization to the case with disorder is well understood. The 
topological quantum number in IQHE, which character- 
izes the quantized Hall conductivity, is known as the first 
Chern number [also refered to as the Thouless-Kohmoto- 
Nightingale-den Nijs (TKNN) integer] and is closely re- 
lated to Berry's phase. In the presence of disorder, the 
TKNN integers defined for a clean system can be general- 



ized. By introducing generalized periodic boundary con- 
ditions and averaging over different boundary condition 
phase, an invariant expression can be constructed appli- 
cable to the situation where many-body interaction and 
substrate disorder are present [TH H3] • Since the bound- 
ary condition phases can be transferred to the Hamil- 
tonian by a unitary transformation, such considerations 
are actually equivalent to thinking of the system as a 
super-cell of an infinite system. Since the infinite system 
which is periodic in the super-cell has translation sym- 
metry, the wave vectors can be well defined, which in 
fact correspond to the boundary condition phases. The 
advantage of such consideration is that since the band 
structure is recovered, one can use the known methods 
to calculate the topological invariant for the infinite sys- 
tem. The finite system under consideration, which is 
now a super-cell of the infinite system, shares the same 
topological properties as the infinite system. This idea 
has been used to study the phase transition in the pres- 
ence of disorder in 2D quantum spin Hall system (QSHE) 
|14) . where the authors found that a metallic region al- 
ways appears between ordinary and topological insulator 
in spin-orbit coupled systems with disorder when there 
is no extra conservation law (Recently the same result is 
also obtained via C*- Algebras [TS]). 

In this paper, we study the topological invariant of the 
disordered 3D system by viewing it as a super-cell of an 
infinite system. We start from a trivial insulator with 
spin-orbit coupling and find that when strong enough 
disorder is introduced into the system, a gap appears at 
half filling and the corresponding topological invariant 
is non-trivial, confirming the existence of the disorder- 
induced non-trivial topological phase. 

To be concrete, we consider a model describing itiner- 
ant electrons with spin-orbit coupling on a cubic lattice 
with the Hamiltonian in the momentum space |161 I17j , 
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FIG. 1: (Color online)The cubic lattice on which the Hamil- 
tonian is defined. As an example, the super-cell has a size 2 3 
(red closed circles). The sites which have connections with 
next (previous) super-cell obtain phases e l ^ i (e - * 1 **), where 
4>i = or 7r and i — x,y,z depends on the TRIM and the 
direction of the bonds. 
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where d±(k) = di(k) ± «d 2 (k), c?o(k) = e — 2i^ i cos£; i , 
cii(k) = — 2Asin(fci) and rf 4 (k) = 27(3 — J2i cosfci) 
(i = 1,2,3). This Hamiltonian is a lattice version of the 
effective low-energy Hamiltonian describing the physics 
of insulators in Bi2Se^, family 18-21J. At half filling, 
depending on the parameters, the system can be a triv- 
ial insulator or a topological insulator. To simulate the 
effects of disorder we consider a random on-site potential 
TijUj^j^j, with Uj uniformly distributed in the range 
(— E/o/2, Uq/2). This kind of disorder respects time- 
reversal symmetry. The lattice on which the Hamiltonian 
resides is a cubic lattice with a finite size L x x L y x L z (the 
lattice constant a = 1). There is no translation symmetry 
in the system when disorder is present. However taking 
the system as a super-cell of an infinite system, transla- 
tion symmetry is recovered. The lattice vector becomes 
a', = Li (the size of the super-cell in i = x,y,z direc- 
tion) and the corresponding reciprocal lattice vectors are 
b 2zl 

1 ~ Li ■ 

In 3D, the topological invariant consists of four Z2 
numbers forming an index (vq; v\Viv^\ which distinguish 
16 topological classes [TOl [22]. Usually it is a difficult 
problem to evaluate them for a given band structure. 
However in the presence of inversion symmetry, the prob- 
lem can be greatly simplified. It has been shown that 
they can be determined from the knowledge of the par- 
ity ^2m(^i) of the 2m-th occupied energy band at the 8 
time reversal invariant momenta (TRIM) I\ that satisfy 
Ti = Ti + G. The 8 TRIM can be expressed in terms 
of primitive reciprocal lattice vectors as ^i=( ni n 2 n 3 ) = 



(nibi + 7i2b2 + n3b3)/2, with rij = 0, 1. Then v a is de- 
termined by the product (—1)"° = rin =o 1 ^n 1 n 2 n 3 , and 
(_ 1 )^=i,2, 3 _ H nj ^ =0il . nt=1 6 ninan3 , where the parity 

product for the occupied bands <J, = Ilm=i &m(Fi)- To 
take advantage of the simplification, we only consider dis- 
order configurations with inversion symmetry. For large 
enough super-cell such consideration will not change the 
underlying physics. So to calculate the topological in- 
variant, we only need to consider the Hamiltonian at the 
eight TRIM and they are equivalent to those of a finite 
system with boundary conditions which are periodic up 
to phases (j> x , <p y , (j) z =0 or tt for boundary sites, as shown 
in Fig. |1| 

To understand what happens when introducing the 
super-cell, we first calculate the Z 2 topological invari- 
ant of a clean (1;000) strong topological insulator. The 
results are shown in Fig. [2] (a) . At half filling, we find 
that 6 = -1 at the T point and S = 1 at other TRIM. 
This result is consistent with the ordinary band structure 
calculations and the reason is explained below. Taking 
a super-cell means enlarging the original unit cell in the 
ai,a2,a3 directions (ai,a2,a3 are the lattice vectors of 
the original lattice). The resulting new Brillouin zone 
(BZ) is folded in the corresponding directions and shrunk 
in size. The occupied states at the eight TRIM for the 
new BZ contain those at the eight TRIM for the origi- 
nal BZ. Though more states which are at other momenta 
of the original BZ will reside on the TRIM of the new 
BZ, they do not change the parity product for the oc- 
cupied states. Thus the product of all eight Ss from the 
super-cell calculation still yield the same 'strong' index 
of the Z 2 topological invariant. However the 'weak' in- 
dex cannot be obtained from this method if the super- 
cell contains an even number of unit cell in the corre- 
sponding direction. As mentioned above, the weak in- 
dices v\ , V2 , vz are the product of four <5s which are in the 
planes ki = 7r,/j 2 = tt, &3 = tt respectively. Suppose that 
the super-cell contains an even number of unit cell in the 
ai direction, then the k\ = and k\ = tt planes will col- 
lapse onto the k[ — plane in the new BZ. Thus the weak 
index determined by the four TRIM on the k[ = n plane 
in the new BZ must be 0. From another point of view, 
this is understandable because the 'weak' indices are re- 
lated to layered 2D quantum spin-Hall states and the 3D 
super-cell naturally fails in calculating the quantities re- 
flecting the 2D physics. When considering systems with 
disorder, the weak indices are eliminated and only the 
strong index remains robust. The topological invariant 
for the system with disorder is the 'strong' index of the 
Z2 topological invariant for clean 3D band insulators. In 
the following, we simply call it the topological invariant 
for the 3D systems with disorder. 

Now we take into account disorder as described above 
and start from a trivial insulator. Fig. [2] (b) shows the 
result of such a calculation on a 8 3 lattice with Uq = 150 
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FIG. 2: (Color online) The parity product for the occupied states at the eight TRIM for a 8 3 super-cell without disorder (a) 
and with disorder (b) . The crosses mark the eigenvalue of the system and the corresponding parity product for the filling up 
to this energy level. Here we use parameters: t — 24meV, A = 20meV, 7 = 16 meV and (a) e = 134meV, corresponding to 
m — —10 meV, where the system is a (1; 000) strong topological insulator; (b)e = 145 meV and Uo = 150 meV, corresponding 
to m = 1 meV, where the clean system is a trivial insulator. The red lines in both figures show the gap range which appears 
at half filling. 
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FIG. 3: (Color online)Distributions of the number of i?2-odd 
band pairs from 500 disorder realizations. The bar heights is 
the fraction of disorder realizations that have a given number 
of band pairs with Z2 = 1 (The numbers are all odd here). 
The parameters are the same as those in Fig. [51(b). 



meV. Generally the disorder will eliminate all degenera- 
cies except those protected by time reversal symmetry. 
So at the TRIM, each eigenvalue of the Hamiltonian 
is doubly degenerate. Though the disorder makes the 
energy spectrum more continuous, a small gap remains 
clear at half filling. Here the gap location is shifted to the 
high energy, which is different from the case in the clean 
system where the gap appears symmetrically at E = 0. 
The parity eigenvalue products S = — 1 at the T point 
and 5 = 1 at other TRIM for half filling show that the 
topological invariant of the system is 1, so the system 
will exhibit non-trivial topological properties. We have 
extended the above calculation to 500 different disorder 
realizations. In all 500 disorder realizations, the topolog- 
ical invariant is 1 for half filling, which further confirms 
the topological phase in the system. Thus through calcu- 
lating the topological invariant of the system, we obtain 
further confirmation of the topological phase induced by 
disorder. This new phase has been termed as 'strong 



topological Anderson insulator' (STAI) [8]. Due to the 
non-trivial topological properties in STAI, topologically 
protected surface states will appear at the surfaces as 
the case in the 'strong topological insulator' phase of the 
clean system, which has already been verified by trans- 
port calculations. 

It is also interesting to look at the number of ^-odd 
band pairs. In 2D QSHE, the number of Z 2 -odd Karmers 
pairs increases linearly with the system size in the metal- 
lic region and grows slowly for the topological insulator 
|14) . Here in 3D, we find a similar result. In Fig. |3j 
we show the distribution of the number of Z2-odd band 
pairs at three lattice sizes 4 3 , 6 3 and 8 3 . We find that 
the location of the mean roughly scales with L 3 (L is the 
lattice size). A detailed analysis on the limited data also 
shows that the number growth is somewhat slower than 
linear. With enough data and larger sizes, finite-size scal- 
ing could be carried out. However this is not accessible 
at present due to the limited computer resources. 

For different disorder realizations, though there still 
exits a gap at half filling, its position on the energy axis 
changes randomly. We extract the energy levels for the 
most energetic electron at half filling and the one just 
above it and show their distributions in Fig. [4j The 
curves in Fig. [4] show Gaussian-like distributions. The 
distributions for half filling and half plus one at some 
TRIM have no overlaps while on the other TRIM have 
overlaps. This means there is no 'true' gap range even 
for the eight TRIM. However the overlaps are already 
very small for the lattice size of 8 3 . We also performed 
the same calculations on lattice sizes of 4 3 and 6 3 . With 
the increasing lattice size the widths of the distributions 
decrease. We therefore attribute the overlaps for the 
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present lattice size to finite size effect and expect them 
to diminish for larger lattice sizes. 

For the parameters used the physical gap size is deter- 
mined by gap at k = (0,0,0) TRIM. We can obtain the 
gap value from the peaks of the distributions and this 
value will approximate the one for the larger lattice size. 
We have determined the phase diagram in the Uq — Ep 
plane from calculations on a 6 3 system and the result is 
shown in Fig. [5] The range of STAI in the Uq — Ep plane 
has also been obtained from conductivity calculations 
and the SCBA, where the weak-disorder boundary fits 
well with each other while the strong-disorder boundary 
does not [5] . It was found that the weak-disorder bound- 
ary marks the crossing of a band edge and the strong- 
disorder boundary marks the crossing of a mobility edge, 
and The SCBA does not work for the strong-disorder 
boundary [TJ. 

The result in Fig. [5] shows that the phase boundaries 
obtained from the present calculation fits well for both 
the weak-disorder and strong-disorder phase boundaries 
determined from the conductivity calculation [8 . The 
range of STAI phase should be determined by the gap 
between the two edges containing extended states. We 
believe that we have extracted such a gap. The reason is 
that the localized states are closely related to the disor- 
der realizations and can be removed by an average over 
disorder realizations. In the phase diagram, increasing 
the strength of disorder, the system first undergoes a 
sharp phase transition, which happens at a critical dis- 
order strength where the gap closes. Then the system 
enters a stable topological phase for a wide range of dis- 
order strength. Finally the system experiences another 
phase transition to a diffusive metal. Close to the latter 
transition, the finite system under consideration is in a 
mixed phase, where some fraction of disorder realizations 
yield a non-trivial topological invariant. 

In conclusion, we have generalized the method of cal- 
culating the topological invariant in disordered 2D QSHE 
to 3D disordered systems. In this method, the finite 3D 
system is viewed as a super-cell of a large lattice with well 
defined wave vector, which allows us to directly use the 
definition of topological invariants for clean band insula- 
tors. The obtained topological invariant can be thought 
of as describing the finite system with disorder. Using the 
inversion symmetry preserving disorder configurations, 
we carried out explicit calculations on a model Hamilto- 
nian (describing the physics of insulators in BiiSe?, fam- 
ily) with on-site disorder. We found that the topological 
invariant for a system that is a trivial insulator in the 
absence of disorder can become nontrivial when strong 
enough disorder is introduced. This result confirms the 
existence of the strong topological Anderson insulator, a 
topological phase in three space dimensions whose exis- 
tence is fundamentally dependent on disorder |S]. 

As additional application of this method we counted 
the number of Z2-odd band pairs. Though this number 




FIG. 4: (Color online) Distributions of energy levels at half 
filling and half plus one filling from 500 disorder realizations. 
The system sizes are 4 3 (blue), 6 3 (green) and 8 3 (red). The pa- 
rameters are the same as those in Fig. |2](b). 
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FIG. 5: (Color online) Phase boundaries for the STAI phase 
in the Uq — Ef plane obtained from the super-cell calcuations. 
The red crosses represent the peaks in energy level distribu- 
tions for half and half plus one fillings and are obtained from 
calculations of 500 disorder realizations on 6 3 system. The 
phase boundaries from the SCBA are also shown for compar- 



varies with the disorder realization, the total number is 
always odd when the system is in the STAI phase and 
it scales roughly linearly with the number of the lattice 
sites. Finally, we obtained the phase boundaries for the 
STAI phase by extracting the gap range, which fits well 
with the known results. 

In the presence of on-site disorder, the diagonal el- 
ements of the Hamiltonian matrix are random, drawn 
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from a statistical distribution. The Hamiltonian can be 
regarded as a random matrix which can be studied in 
the framework of the random matrix theory (RMT). The 
RMT has already been applied to IQHE and disordered 
superconductors 23-25J. Wc expect some insights into 
the present problem from RMT, which we leave to future 
work. 

Acknowledgment. — Authors are indebted to I. Garate, 
M. Franz, G. Refael, G. Rosenberg, C. Weeks and M. 
Vazifeh for stimulating discussions. Support for this work 
came from NSERC, ClfAR and The China Scholarship 
Council. 



[1] B.A. Bernevig, T.L. Hughes, and S.-C. Zhang, Science 

314 1757 (2006). 
[2] J.E. Moore, Nature 464, 194 (2010). 
[3] M.Z. Hasan, C.L. Kane, arXiv: 1002.3895. 
[4] X.-L. Qi, T.L. Hughes, and S.-C. Zhang, Phys. Rev. B 

78 195424 (2008). 
[5] J. Li, R.-L. Chu, J. K. Jain, and S.-Q. Shen, Phys. Rev. 

Lett. 102 136806 (2009). 
[6] H. Jiang, L. Wang, Q.-F. Sun, and X.-C. Xie, Phys. Rev. 

B 80 165316 (2009). 
[7] CW. Groth, M. Wimmer, A.R. Akhmerov, 

J. Tworzydlo, and C.W.J. Beenakker. Phys. Rev. 

Lett. 103 196805 (2009). 
[8] H.-M. Guo, G. Rosenberg, G. Refael, and M. Franz, 

larXiv:1006.2777l 
[9] J. E. Moore and L. Balents, Phys. Rev. B 75 121306(R) 



(2007). 

[10] L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007). 
[11] D.J. Thouless, M. Kohmoto, M.P. Nightingale, and M. 

den Nijs, Phys. Rev. Lett. 49 405 (1982). 
[12] Q. Niu, D.J. Thouless, and Y.-S. Wu, Phys. Rev. B 31, 

3372 (1985). 

[13] J.E. Avron and R. Seiler, Phys. Rev. Lett. 54, 259 (1985). 
[14] A.M. Essin and J.E. Moore, Phys. Rev. B 76, 165307 
(2007). 

[15] T.A. Loring and M.B. Hastings, arXiv: 1005.4883. 
[16] P. Hosur, S. Ryu, and A. Vishwanath, Phys. Rev. B 81, 
045120 (2010). 

[17] G. Rosenberg and M. Franz, Phys. Rev. B 82, 035105 
(2010). 

[18] H. Zhang, C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.- 

C. Zhang, Nature Phys. 5, 438 (2009). 
[19] Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, 

A. Bansil, D. Grauer, Y. S. Hor, R. J. Cava and M. 

Z. Hasan Nature Phys. 5, 398 (2009). 
[20] Y. L. Chen, J. G. Analytis, J.-H. Chu, Z. K. Liu, S.- 

K. Mo, X. L. Qi, H. J. Zhang, D. H. Lu, X. Dai, Z. Fang, 

S. C. Zhang, I. R. Fisher, Z. Hussain, and Z.-X. Shen 

Science 325, 178-181 (2009). 
[21] D. Hsieh, Y. Xia, D. Qian, L. Wray, F. Meier, J. H. Dil, 

J. Osterwalder, L. Patthey, A. V. Fedorov, H. Lin, 

A. Bansil, D. Grauer, Y. S. Hor, R. J. Cava and M. 

Z. Hasan, Phys. Rev. Lett. 103, 146401 (2009). 
[22] H.-M. Guo and M. Franz Phys. Rev. Lett. 103, 206805 

(2009). 

[23] David J. E. Callaway Phys. Rev. B 43, 8641 (1991). 
[24] Martin Janssen and Krystian Pracz, Phys. Rev. E 61, 
6278 (2000). 

[25] Sail R. Bahcall, Phys. Rev. Lett. 77, 5276 (1996). 



